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i ' Abstract 

Tomograms, a generalization of the Radon transform to arbitrary 

pairs of non-commuting operators, are positive bilinear transforms 

cZ2 ' with a rigorous probabilistic interpretation which provide a full char- 

. ,— i ■ acterization of the signal and are robust in the presence of noise. We 

^. provide an explicit construction of tomogram transforms for many 

JS^ \ pairs of noncommuting operators in one and two dimensions and il- 

i ^~ *' lustrations of their use for denoising, detection of small signals and 

component separation. 

> 

(N ! 1 Introduction 

Integral transforms [T] [2] are very useful for signal processing in communi- 
cations, engineering, medicine, physics, etc. Linear and bilinear transforms 
have been used. Among the linear transforms, Fourier [3] and wavelets jl] 
[5] j6] are the most popular. Among the bilinear ones, the Wigner-Ville qua- 
sidistribution [7] [8] provides information in the joint time-frequency domain 
/s^ with good energy resolution. A joint time-frequency description of signals 

is important, because in many applications (biomedical, seismic, radar, etc.) 
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the signals are of finite (sometimes very short) duration. However, the os- 
cillating cross-terms in the Wigner-Ville quasidistribution make the inter- 
pretation of this transform a difficult matter. Even if the average of the 
cross-terms is small, their amplitude may be greater than the signal in time- 
frequency regions that carry no physical information. To profit from the 
time-frequency energy resolution of the bilinear transforms while control- 
ling the cross-terms problem, modifications to the Wigner-Ville transform 
have been proposed. Transforms in the Cohen class [9] [10] make a two- 
dimensional filtering of the Wigner-Ville quasidistribution and the Gabor 
spectrogram [11] is a truncated version of this quasidistribution. The diffi- 
culties with the physical interpretation of quasidistributions arise from the 
fact that time and frequency correspond to two noncommutative operators. 
Hence a joint probability density cannot be defined. Even in the case of 
positive quasiprobabilities like the Husimi-Kano function [12] [13], an inter- 
pretation as a joint probability distribution is also not possible because the 
two arguments of the function are not simultaneously measurable random 
variables. 

Recently, a new type of strictly positive bilinear transforms has been 
proposed [H] [15], called tomograms, which are a generalization of the Radon 
transform [16] to noncommutative pairs of operators. The Radon- Wigner 
transform [T7] [18] is a particular case of such noncommutative tomography 
technique. The tomograms are strictly positive probability densities, provide 
a full characterization of the signal and are robust in the presence of noise. 

A unified framework to characterize linear transforms, quasidistributions 
and tomograms was developed in Ref.[15]. This is briefly summarized in 
Section 2. Then Sections 3,4,6 and 7 contains an explicit construction of 
tomogram transforms for many pairs of noncommuting operators in one and 
two dimensions. Some of these transforms have been used in the past [22] 
[27], others are completely new. 

It is in the time-frequency plane that most signal processing experts have 
developed their intuition, not in the eigenspaces associated to the new to- 
mograms. Therefore, to provide a qualitative intuition on the way the tomo- 
grams explore the time-frequency plane, we have provided graphical spectro- 
grams of the eigenstates on which the signal is projected by the tomograms. 

In Section 5, an interpretation of the tomograms is given as operator sym- 
bols of the set of projection operators in the space of signals. This provides 
a very general framework to deal with all kinds of custom-designed integral 
transforms both for deterministic and random signals. It also provides an 



alternative framework for an algebraic formulation of signal processing. 

Finally, an illustration of how such transforms may be used to analyze 
signals is contained in Section 8. A brief review of denoising, detection of 
small signals and component separation, done in the past, is included as well 
as an application of one of the new transforms. 

2 Linear transforms, quasi-distributions and 
tomograms 

Consider signals f(t) as vectors | /) in a dense nuclear subspace M of 
a Hilbert space 1-L with dual space M* (with the canonical identification 
A" C A/"*) and a family of operators {U(a) : a G I, I C M n } defined on M* . 
In most cases of interest U (a) generates a unitary group U (a) = e tB ^ a \ In 
this setting three types of integral transforms are constructed. 

Let h G A/"* be a reference vector and let U be such that the linear span 
of {U(a)h G Af* : a G /} is dense in M* . In the set {U(a)h}, a complete set 
of vectors can be chosen to serve as a basis. 

1 - Linear transforms 

wf\a) = (U{a)h\f) (1) 

2 - Quasi-distributions 

Q f (a) = (U(a)f\f) (2) 

3 - Tomograms 

If U (a) is a unitary operator there is a self-adjoint operator B (a) such 
that U (a) = e lB ^ . The tomogram is 

Mf\x) = {f\5{B{a)-X)\f) (3) 

X takes values on the spectrum of B(a). Considering a set of generalized 
eigenstates (in A/"*) of B(a), one obtains for the kernel 

(Y | 5 (B(a) - X) | Y') = S(Y' - X) S(Y - Y') = (Y\X)(X\ Y') 

Therefore, we may identify 5 {B(a) — X) with the projector | X)(X | 

5(B(a)-X)=\X)(X\=P x 



From this, it follows 

M\ B) = (f \8(B(a)-X) | f) = (f | X)(X | /) = \(X \ f)\ 2 (4) 

showing the positivity of the tomogram and its nature as the squared ampli- 
tude of the projection on generalized eigenvectors of B(a). Let, by a unitary 
transformation S, B(a) be transformed to 



SB(a)£fl = B'\ 



a 



If {| Z)} is the set of (generalized) eigenvectors of B (a), {S^ \ Z)} is a set 
of eigenvectors for B. Therefore, 

Mf\z) = (f | S(B(a) - Z) | /) = \{Z | S \ f)\ 2 = (f \ S^ \ Z)(Z \ S \ f) 

For normalized | /), 



(f\f) 



the tomogram is normalized 



J M { f B) (X) dX = 1 (5) 



It is a probability distribution for the random variable X corresponding to the 
observable defined by the operator B (a). The tomogram is a homogeneous 
function 

Mf tv) [X) = \p\Mf\pX) (6) 

Examples: 

If U (a) is unitary generated by Bp (a) = ait + ia.2J[ t and h is a (gen- 
eralized) eigenvector of the time-translation operator the linear transform 
Wf (a) is the Fourier transform. For the same Bp (a), the quasi-distribution 
Qf(a) is the ambiguity function. 

The Wigner-Ville transform [7] [8] is the quasi-distribution Qf(a) for the 
following B— operator 

(j2 _ d?_ _ j 
B^ wv \ ai , a 2 ) = -i2 ai ± - 2a 2 t + -A f '- (7) 

The wavelet transform is W f (a) for Byy (a) = a\D + ict2-^, D being 
the dilation operator D = — | (it4: + i4rt). The wavelets h S)T (t) are kernel 



functions generated from a basic wavelet h(r) by means of a translation and 
a rescaling (—00 < r < 00, s > 0): 



, / n 1 , /* — r 



(8) 



using the operator 

£/"' '(r, s) = exp(zra)) exp(ilog sD), (9) 

M*) = tf (A)t fa*)M*)- (10) 

For normalized h(t) the wavelets h SyT (t) satisfy the normalization condition 

\h S:T (t)\ 2 dt = 1. 

The basic wavelet (reference vector) may have different forms, for example, 

1 



h(t) 



e ia, t e -t 2 /2 
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or 



(11) 
(12) 



/ l (t) = (l-t 2 )e- t / 2 

called the Mexican hat wavelet. 

The Bertrand transform [19] [20] is Qf{ot) for 5^- 

Linear, bilinear and tomogram transforms are related to one another by 



/f >, 



W'(X) = T JQ, 



(fcfl), 



a) e 



-ifcV 



r/A' 



Qf\a)= Mf /p \x)e ipX dX 
Qf\a) = W { /\a) 



Wf\a) = -\e 



AX 



-M { *\X) + IM { *\X) 



dX 



with 



h) =1 h)+ I /}; I f 3 ) =| h)- I /> 
/ a )=IM + »l/>; \h)=\h)-i\f) 



3 One- dimensional tomograms 

As shown in (J3J) a tomogram corresponds to projections on the eigenstates 
of the B operators. These operators are linear combinations of different 
(commuting or noncommuting) operators, 

B = ^O x + v0 2 

Therefore the tomogram explores the signal along lines in the plane (Oi, O2). 
For example the tomogram 

Mf ] (X, f i,v) = (f\6 (fit + uu-X)\f) (13) 

with ui = i4, is the expectation value of an operator delta- function in the 
state I /), the support of the delta-function being a line in the time-frequency 
plane 

X = fit + vu (14) 

Therefore, M\ (X, n,u) is the marginal distribution of the variable X along 
this line in the time-frequency plane. The line is rotated and rescaled when 
one changes the parameters /x and v. In this way, the whole time-frequency 
plane is sampled and the tomographic transform contains all the information 
on the signal. 

It is clear that, instead of marginals collected along straight lines on the 
time-frequency plane, one may use other curves to sample this space. It has 
been shown in [15] that the tomograms associated to the affine group, for 
example 

M^ (X, fi, u) = (f I 5 Lt + v l ^± - x) I /) (15) 

correspond to hyperbolas in the time-frequency plane. This point of view has 
been further explored in [21] defining tomograms in terms of marginals over 
surfaces generated by deformations of families of hyperplanes or quadrics. 
However not all tomograms may be defined as marginals on lines in the 
time-frequency plane. 

Here we construct the tomograms corresponding to a large set of operators 
defined in terms of (one-dimensional) time. Of particular interest are the 
tomograms associated to finite-dimensional Lie algebras. 



3.1 ID conformal group tomograms 

The generators of the one- dimensional conformal group are, 

<" = <* 

D = ihi + i) (16) 

One may construct tomograms using the following operators: 
Time-frequency 

B 1= i2t + iiy— (17) 



Time-scale 
Frequency-scale 
Time- conformal 



(It 



B 2 = ,,t + I u(tj t + ^) (18) 



d ( d 1 , 

B ^ i "it + iv { t Jt + ^ (19) 



B A = /it + iv(t 2 — + t) (20) 



The construction of the tomograms reduces to the calculation of the gen- 
eralized eigenvectors of each one of the Bi operators 
B x fa (//, v, t, X) = Xfa (//, is, t, X) 

( \xt 2 tX\ . . 

fi(H,v,t,X) = expi[- (21) 

\2is is J 

with normalization 

dtfa (/i, v, t, X) fa fa, is, t, X') = 2tvis5 (X - X') (22) 



B 2 fa (//, v, t, X) = Xfa (fi, v, t, X) 



fa (a*, v, t, X) = — = exp i ( log \t\ ) (23) 

\t\ \is is 



I dti)* 2 (/i, is, t, X) fa (/i, is, t, X') = Attis5 (X - X') (24) 



B 3 ip 3 fa, v, w, X) = X^ 3 (/i, i/, w, X) 



ip 3 (/i, i/, t, X) = exp (-i) i-u log \u\ 

V V 



dunpl (/i, v, co, X) tpi (/i, v, cu, X') = 2ttu5 (X — X') 
£4^4 (//, v, t, X) = X^ 4 (/i, ^, t, X) 



1 ■ f x V 1 , 

t/> 4 (/i, z/, t, X) = — expz — + -log £ 

|r| \z/£ ^ 



dtipl (/i, i/, t, s) ?/> 4 (/i, ^, t, s') = 2nv5 (s — s') 



Then the tomograms are: 
Time-frequency tomogram 



2ttM 



M 1 (ii,v,X) = 

Time-scale tomogram 

M 2 (fi,u,X) 
Frequency-scale tomogram 
M 3 (fi,u,X) 



cxp 



i\xt 2 itX 
~2v V 



f(t)dt 



2n\u\ 



dt 



M P K^-fi°gi*i)] 



2ir\u\ 



/M e H(£u,-£ logH)] 



gL> ; _ e 



/(w) being the Fourier transform of /(£) 
Time-conformal tomogram 



M 4 (/i,z/,X) 



1 



2ttM 



M P Ki+fio g i*i)] 



dt -TT-e 



(25) 
(26) 

(27) 
(28) 



(29) 



(30) 



(31) 



(32) 



The tomograms Mi, M 2 and M 4 interpolate between the (squared) time signal 
and its projection on the ipt (/1, v, t, X) functions for /i = 0. FigJT] shows the 
typical behavior of the real part of these functions. 
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Figure 1: Typical behaviour of the real part of the functions 1^1,1^2 and ^4 
at /i = 



Figsj2j3] and H] illustrate how the tomograms M\,M2 and M4 explore the 
time-frequency space by plotting the spectrograms of typical vectors ipi,ip2 
and -04. 

In a similar way, tomograms may be constructed for any operator of the 
general type 

. / . . d ldgCt)' 



the generalized eigenvectors being 

if) g (ji,u,t,X) = |#(£)r 1/2 expi 



X 



ds u 

+ 



g(s) v] g(s) 



sds 



3.2 Another finite-dimensional algebra 

Another finite-dimensional Lie algebra which may be used to construct to- 
mograms, exploring other features of the signals, is generated by 1, t and 
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Figure 2: Modulus of the Short-time Fourier transform of 4 vectors of the 
time-frequency tomogram for some fixed 9, \x — cos 9, v = sin9. A vector 
is a linear chirp, hence a line in the time-frequency plane. Moreover, each 
vector is a frequency-translated version of the one which starts at the origin. 
Since it forms an orthogonal basis, the sum of all the vectors cover the entire 
time-frequency plane. The parameter 9 allows to change the slope of the line 
in the time-frequency plane. 
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Figure 3: Modulus of the Short-time Fourier transform of 4 vectors of the 
time-scale tomogram for /j, — 0, v = 1 (left) and fj, = y(2)/2, v = a/(2)/2 
(right). Each vector is an hyperbolic chirp. Two of them correspond to 
positive X and two of them to negative X. Due to the sampling used in the 
numerical computation, some aliasing phenomenum occurs at times close to 
zero. There is a axis of symetry: the line of zero frequency on the left graph. 
This axis is shifted in frequency when \i and v are changed. 



ui - 
D 



a 



% dt 



H + \) 



t 2 + i 



i / d 2 

2 \dP 



Of special interest are the tomograms related to the operators 

B F = jit + vF 



and 



B a = fit + z/cr 



As before, the construction of the tomograms relies on finding a complete 
set of generalized eigenvectors for the operators Bp and B a . With y = t + - 
one defines creation and annihilation operators 



1 

72 
1 

71 



y + 
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dy 

d 
dy 
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Figure 4: Modulus of the Short-time Fourier transform of 4 vectors of the 
time-conformal tomogram for fi — 0, v — 1. Due to the sampling used in 
the numerical computation, some aliasing phenomenum occurs at times close 
to zero. Some interferences between the vectors occur for large time. Two 
vectors correspond to positive X and two to negative X. 
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B„ = v I aa 



obtaining 

2z/ 2 

2z/ 2 
Therefore for Bp one has an orthonormalized complete set of eigenvectors 

C } (*) =u n (t+ £) 
with a discrete set of eigenvalues X n = v (n + |) — ^ 

5^ (*) = ^n^i F) (t) 

the function u n being 

u n (y) = (n^n\y 1/2 (y - £)" ' e~4 
The tomogram M j- (//, u, X n ) is 



r(^) 



My>{n,v,X n 



ri F > (t) f (t) dt 



My) 



For B a one uses a basis of coherent states 

„Aa' — A*a 
o 

wi^ A" 

;i 

with decomposition of identity 



«o(y) 

t^ vn! 



1 



7T 



A (l/)^(l/)rf 2 A = l 



Then, a set of generalized eigenstates of B a is 

^f (/i, M) = A (t + ^) 

with eigenvalues 

B a ^ ) {n,v,t)=X x ^ ) { l i,v,t) 
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Xx = u[\ 



JL. 

2z/ 2 



the tomogram being 



Mfi^v.X, 



(o)* 



W(/M,V,t)f{t)dt 



This tomogram is closely related to the Sudarshan-Glauber P-representation 

[23 |23]. 

4 Multidimensional tomograms 

Several types of multidimensional tomograms may be obtained from gen- 
eralizations of the one-dimensional ones. Consider a signal f{ti,t 2 ). The 
tomogram will depend on a vector variable X = (Xi,X 2 ) and four real pa- 
rameters /ii, ^2i v ii and v 2 . For example, the two-dimensional time-frequency 
tomogram will be 



M(X,ft,u) 



1 



1 



/(ti,t 2 )exp ( % -^t\ - % ^t x + ^t 2 2 -^t 2 ) dhdh 
2z/i v\ 2u 2 v 2 



2 r K\Vx\ 2ll\l> 2 

From this one may also construct a center of mass tomogram 



(33) 



M cm {Y, fl,v)= / M(X,jl, V) 5{Y - X, - X 2 ) dX, dX 2 = 6{Y-X 1 - X, 



1 



2ti\vi\ 2ix\v 2 \ 



x 



. l/Xl 2 Z«iXi Z/i 2+2 122-^2 

f(ti,t 2 )dt 1 dt 2 exp ( - — *! + 7^ — 1 2 

2V\ V\ lv 2 v 2 



dX 1 dX 2 



the center of mass tomogram being normalized 



f M cra (X,fI,u)dX = l 



and a homogeneous function 



1 



M cm (XX, A/Z, Xu) = — M cm (X, //, z7) 
II 

The generalization to N channels is straightforward. 
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As in the one- dimensional case, useful tomograms may be constructed 
from the operators of Lie algebras. For example, given the generators of the 
conformal algebra in IR d , d > 2, 



u k = i m 



H, 
D = i (t • V + |) 

Rj,k = i \tjai~ ~ ik atJ 

Let, in two dimensions, t\=t and t 2 = x. The tomograms corresponding to 
the operators 

B u = Hit + /jL 2 x + v-lUi + z/ 2 o;2 

Bd = Hit + \i%x + uD 

Buj = Hit + H2X + v\K x + v 2 K 2 

are, as in (|33l) . straightforward generalizations of the corresponding one- 
dimensional ones. For the operator 

Br = Hit + H2X + vR.1,2 
the eigenstates are 

ijr R ' f fj,, v, x, t, X J = exp — I [i\x — H2t + X tan -1 — ] 

and the tomogram 

M f (jl, v,X\= [ i) {R) * (h, v, x, t, x) f (x, t) dxdt 

5 The tomograms as operator symbols 

Tomograms may be described not only as amplitudes of projections on a 
complete basis of eigenvectors of a family of operators, but also as opera- 
tor symbols. That is, as a map of operators to a space of functions where 
the operators non-commutativity is replaced by a modification of the usual 
product to a star-product. 

15 



Let A be an operator in Hilbert space % and U(x), D(x) two families of 
operators called dequantizers and quantizers, respectively, such that 

Trl,U(x)D(x)\ =5(x-x) (34) 

The labels x (with components xi, x%, . . . x n ) are coordinates in a linear space 
V where the functions (operator symbols) are defined. Some of the coordi- 
nates may take discrete values, then the delta function in (l34"j) should be 
understood as a Kronecker delta. Provided the property fl34|) is satisfied, one 
defines the symbol of the operator A by the formula 

f A (x)=Tr{u(x)A}, (35) 

assuming the trace to exist. In view of (1341 . one has the reconstruction 
formula 

A = / f A {x)D(x) dx (36) 

The role of quantizers and dequantizers may be exchanged. Then 



f%{x) = Tr \D{x) A) (37) 

is called the dual symbol of Ja^x) and the reconstruction formula is 

A= f f d A (x)U{x)dx (38) 



Symbols of operators can be multiplied using the star-product kernel as fol- 
lows 



Ja(x) * f B (x) = / f A (y)f B (z)K(y, z, x) dydz (39) 

the kernel being 

K(y,z,x) = Tr {b{y)D{z)U{x)} (40) 

The star-product is associative, 

{f A (x) * f B {x)) * fcix) = /a(x) * (f B (x) * fc(x)) (41) 

this property corresponding to the associativity of the product of operators 
in Hilbert space. 

16 



With the dual symbols the trace of an operator may be written in integral 
form 

Tr { ifi} = J fi(x)f B (x) dx = J f d B (x)f A (x) dx. (42) 

For two different symbols Ja{x) and /^(y) corresponding, respectively, to 
the pairs (U(x),D(x)) and (Ui(y),Di(y)), one has the relation 



Ja(x) = j f A (y)K(x,y)dy, (43) 

with intertwining kernel 

K(x,y) = Tr[b 1 (i/)U(g)} (44) 

Let now each signal / (t) be identified with the projection operator 11^ 
on the function / (t), denoted by 

n/ = |/>(/| (45) 

Then the tomograms and also other transforms are symbols of the projection 
operators for several choices of quantizers and dequantizers. 

Some examples: 

# The Wigner-Ville function: is the symbol of | /)(/ | corresponding to 
the dequantizer 

U(x) = 2V(2a)P, a = l^L (46) 

v2 

where P is the inversion operator 

Pf(t) = f(-t) (47) 

and / D('j) is a "displacement" operator 



15(7) = exp 
The quantizer operator is 



1 a d ^ i * a d 



(48) 



D(x):=D(t,co) = ^-U(t,u), (49) 



t and u being time and frequency. 
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The Wigner-Ville function is 



W(t,u)=2Tt\\f)(f\D(2a)D 



or, in integral form 



W{t,u) = 2 J f*(t)V(2a)f(-t)dt 



(50) 



(51) 



# The symplectic tomogram or time-frequency tomogram of | /)(/ | cor- 
responds to the dequantizer 



with 



U(x) := U(X, fi,v) = 5 (XI - fit - uu) 



tf(t)=tf(t), w/(t) = -.^ /(*) 



and X, fi, v e R. The quantizer of the symplectic tomogram is 
D(x) := D(X, fi, v) = — exp \i (XI — yd, — v0n\ 

# The optical tomogram is the same as above for the case 

ji = cos 9, v = sin 9. 
Thus the optical tomogram is 
M(X,9) -- 



(52) 
(53) 

(54) 

(55) 



Tr {|/)(/ 


\s(xi- 


{d — vCj\ } 








i-|r(t)e^exp 


/ d \ 

ik X — t cos 9 + i — - sin 9 
V dt J 


f(t)dtdk 


1 


ff(t)exp 


". /cot0 2 XtY 
[ l { 2 l sm9j_ 


dt 


2 

(56) 


2tt sin 9 



One important feature of the formulation of tomograms as operator sym- 
bols is that one may work with deterministic signals / (t) as easily as with 
probabilistic ones. In this latter case the projector in ( |4"5]) would be replaced 
by 



n, 



P»\U){U\d/j, 



(57) 
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with J p^dfj, = 1, the tomogram being the symbol of this new operator. 

This also provides a framework for an algebraic formulation of signal 
processing, perhaps more general than the one described in [24] [25]. There, 
a signal model is a triple (A,Ai, $) A being an algebra of linear filters, M. a 
^4-module and $ a map from the vector space of signals to the module. With 
the operator symbol interpretation both (deterministic or random) signals 
and (linear or nonlinear) transformations on signals are operators. By the 
application of the dequantizer (Eq. [35]) they are mapped onto functions, the 
filter operations becoming star-products. 

6 Rotated-time tomography 

Now we consider a version of tomography where a discrete random vari- 
able is used as an argument of the probability distribution function. We 
call this tomography rotated time tomography. It is a variant of the spin- 
tomographic approach for the description of discrete spin states in quantum 
mechanics. For a finite duration signal /(£), with < t < T,we con- 
sider discrete values of time f(t m ) = f m , where with the labeling m = 
—j, — j + 1, — j + 2, . . . , 0, 1, . . . , j — 1, j they are like the components of a 
spinor | /). This means that we split the interval [0,T] onto N parts at time 
values t_j, t_ J+ i, . . . ,tj and replace the signal f(t), a function of continuous 
time, by a discrete set of values organized as a spinor. By dividing by a 
factor we normalize the spinor, i.e., 

(/ 1 /) = E i^i 2 = x ( 58 ) 

m=-j 

Without loss of generality, we consider the "spin" values to be integers, i.e., 
j = 0, 1, 2, . . . and use an odd number N = 2j + 1 of values. 

In this setting, | /) being a column vector, we construct the NxN matrix 

P =\f)(f I (59) 

with matrix elements 

Pmm' = fmfm'- (60) 

The tomogram is defined as the probability-distribution function 

M(m,u) = \(m | u | /)| 2 , m = -j, ...,j - l,j (61) 
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where u is the unitary NxN matrix 

uu ] = ljv (62) 

For this matrix we use an unitary irreducible representation of the rotation 
group (or SU(2)) with matrix elements 



Umm'(9) 



-IV- 



(m + m! 



(j + m)l(j + m') 
(j — m)\(j — m')\ J V 2 
6 



1/2 / /) \ m—m' / /i\ m+m' 



sin - COS 



xJ r j _ m [2m + l,m + m' 2 -) (63) 



J-j-m being a function with Jacobi polynomial structure expressed in terms 
of hypergeometric function as 

F n (a,b,t) = F(-n, a + n,b;t)= ^ ~ * )! ^(1 - t) fe ' a (j^j " [^^(l - t) Q - fe+1 ] 

(64) 
The dequantizer in the rotated-time tomography is 

U(x) = U(m, n) = 8{ml — u*J g u) = 5 ( ml — nJ J (65) 

where J z is the matrix with diagonal matrix elements 

The vector n = (sin 6 cos ip, sin 9 sin 93, cos #) determines a direction in 3D 
space. The matrix (163]) was written for ip = but, if this angle is nonzero, 
the matrix element has to be multiplied by the phase factor e imip . 

The quantizer can take several forms: 

In integral form, it reads 

D(m,n) = / sin 2 — exp(— iJn)^d^{- ■ ■ ) (67) 

The tomogram Ai (m, u) is a nonnegative normalized probability distribution 
depending on the direction n, i.e., M.(m, u) > and 

j 
Y^ M(m,u) = l (68) 

m=-j 
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To compute the tomogram for a given direction with angles <p = and 9, one 
has to estimate 

3 

M(m,9) = Yl <m'(0)fmf™»u m »m(0) (69) 

m" ,m'=—j 

where the matrix u m » m (8) is given by f)63p . The following form for the matrix 
u m i m {9) is more convenient for numerical calculations: 



Um'm\y) 



{j + m')\{j - m')\] 1/2 f 0\ m+m / 0' „,„■-,„.,„■+,„. ,,, 

— — - — cos- sin - P. , ' cos# 

(j + m)\(j-m)\ J V 2/ V 2/ J ~ m 



(70) 
where P°' b are Jacobi polynomials. 

In principle, one could use not only the unitary matrix in ( 1631) but arbi- 
trary unitary matrices. They contain a larger number of parameters (equal 
to iV 2 — 1) and can provide additional information on the signal structure. 

How the time-rotated tomogram explores the time-frequency plane is, as 
before, illustrated by spectrograms of the eigenstates (Figsj5] and [6]). For 
m = 0, formula (70) reduces to the set of normalized associated Legendre 
functions L™ : 



The normalized associated Legendre functions are related to the unmormal- 
ized ones P™ through: 



I2 j + l(j-m') 
2 (j + mf) 



Lf (oosW) = .p^^ii^oos, 



In the tomogram, 9 is the parameter labelling the vectors of the basis as- 
sociated to m = 0,m'. The index j is the variable. In order to illustrate 
the effect of this tomogram, we computed numerically some vectors in the 
time-frequency plane (Figs. |5]and|6]). In the discrete setting, If we choose 
mf = N, where iV is the number of points, the {L^}j form an orthonor- 
mal basis of the discrete time-frequency plane. Hence the projection on the 
eigenvectors of the rotated tomogram with m = 0, mf = N can be seen as the 
projection on the bended lines in the time-frequency plane. This tomogram 
should be adapted for the study of functions which possess certain symetry 
in the time-frequency plane. 
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L5(X) 



i 



1-5 

-U) 

-15 

1-20 



Figure 5: Modulus of the short-time Fourier transform of the sum of 4 
vectors of the rotated-time tomogram. Each u m >^(9,j) is a bended line 
in the time-frequency plane where 9 fix the size. Here N = 1571 and 
9 = 7r/8, 2tt/8, 37r/8, tt/2. As 9 increases, the line is stretched to the right 
until it breaks in two parts for tt/2. 
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Figure 6: Modulus of the short-time Fourier transform of the sum of 4 vectors 
of the rotated-time tomogram. Here N = 1571 and 6 = 57r/8, 67r/8, 7tt/8 and 
Ti — 7r/16. Each u m >fi(9,j) is made of two bended lines in the time-frequency 
plane, one in the upper-half plane and one in the lower-half plane. 
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7 Hermite basis tomography 

Here we consider a dequantizer 

U(n,a) = V(a) | n)(n \ V j (a), a = \a\e Wa 
and a quantizer 



D(n, a) 



IX 



1-A 2 ) V A - 1 



Via) - V -a) 



A + l 



(71) 



(72) 



where — 1 < A < 1 is an arbitrary parameter and n is related to the order 
of an Hermite polynomial. This is analogous to the use of a photon number 
basis in quantum optics. 

For any signal /(£), one has the probability distribution (tomogram) 



M f (n,a)=Tr \f)(f\U(n,a) 
and, from the tomogram, the signal is reconstructed by 

oo „ 

1/X/hE d 2 aM(n,a)D(n,\) 

?1=0 J 

One has Ai(n, a) > and 

J2M f (n,a) = l 

n=0 

for any complex a. For an arbitrary operator A, one has 

oo „ 

IA = J2 d 2 aD(n,a)Tr (u(n,a)A 



n=0 ' 



where I is the identity operator. 

The explicit form of the tomogram for a signal function f(t) is 



M f (n,X)= (f\V(a)\n) 



F(t)f n ,a(t)dt 



where 



f n>a (t) = V(a) Irr^^nir^e-^^H^t) 



(73) 



(74) 



(75) 



(76) 



(77) 



(78) 
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H n (t) being an Hermite polynomial. 
Thus, one has 

/„,«(*) = Tx' ll \2 n n\)- l / 2 e-^- a ^' A e^- a ^ V ~ 2 e- p ' 2 H n it) (79) 
and 

* = «-- ^-- (80) 

For fixed |a| the tomogram is a function of the discrete set n = 0, 1, . . . and 
the phase factor 9 a . 

How the Hermite basis tomogram explores the time-frequency plane is, 
as before, illustrated by spectrograms of the eigenstates (FigJTJ). In the par- 
ticular case where a = 0, the functions f n $ are the Hermite functions. Their 
time-frequency representation has been calculated on Figure It shows 
that the tomogram at a = is suited for rotation invariant functions in the 
time-frequency plane. One can see from (79) that: for real a this pattern 
is shifted in time and for purely imaginary a the pattern is shifted in fre- 
quency. The pattern can be shifted in both time and frequency by choosing 
the appropriate complex value for a. 

8 Some applications 

The tomograms are squared amplitudes of the signal projections on families of 
unitarily equivalent basis (labelled by the /i, v parameters). By inspecting the 
unfolding of these (probability) amplitudes as the parameters change, several 
features of the signals are put into evidence. Here we review briefly three such 
applications, namely denoising, detection of small signals and component 
decomposition, which use the time-frequency tomogram. Then the time-scale 
tomogram will be used to analyse a turbulent velocity fluctuations signal. 

For the finite-time signals, instead of (12 9p . we consider the finite-time 
tomogram 



Mi(0,X) 

with 



to+T 



X 1 ) 



nwiwdt 



to 



2 



2 



</,# 1} > (81) 



1 1 j . 1 ( i cos 9 2 iX 



ip^l (t) = —= exp 1 2 1 (82) 

re,x\J ^ v \2€mQ sin# ' 
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SOU 1DDU J 500 

I imir (..->amplcsl 

Figure 7: Modulus of the short-time Fourier transform of the sum of 4 Her- 
mite functions. Each ring is a Hermite function. Here, the number of points 
is N = 2000. The picture has been centered, the origin has been set to Time 
t = 1000, Frequency / = 0. That is to say, t = -N/2 + I At for I E [0, N), 
At = 1. The smallest circle is for n = 5 and in increasing size order n = 500, 
n = 1000, n = 1500, respectively. 
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and (j, = cos 9,u = sin 9. 

9 is a parameter that interpolates between the time and the frequency 
operators, running from to n/2 whereas X is allowed to be any real number. 

8.1 Detection of small signals 

As an example [14] consider a signal generated as a superposition of a nor- 
mally distributed random amplitude - random phase noise (with total dura- 
tion T = 1) with a sinusoidal signal of same average amplitude, operating 
only during the time 0.45 — 0.55. The signal to noise power ratio is 1/10. The 
true nature of the signal is not revealed neither from its time development 
nor from its Fourier spectrum. However computing the tomogram (see the 
contour plot in FigJH]) one sees clearly a sequence of small peaks connecting 




Figure 8: Detection of small signals in noise 
a time around 0.5 to a frequency around 200. The signature that the signal 
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leaves on the tomogram is a manifestation of the fact that, despite its low 
signal to noise ratio, there is a certain number of directions in the (t, u) plane 
along which detection happens to be more favorable. For different trials the 
coherent peaks appear at different locations, but the overall geometry of the 
ridge is the same. On the other hand, a ridge of small peaks is reliable be- 
cause the rigorous probability interpretation of M(8, X) renders the method 
immune to spurious effects. 

8.2 Denoising and component decomposition 

Most natural and man-made signals are nonstationary and have a multicom- 
ponent structure. Therefore separation of its components is an issue of great 
technological relevance. However, the concept of signal component is not 
uniquely defined. The notion of component depends as much on the observer 
as on the observed object. When we speak about a component of a signal 
we are in fact referring to a particular feature of the signal that we want to 
emphasize. For signals that have distinct features both in the time and the 
frequency domain, the time-frequency tomogram is an appropriate tool. 

Here again consider finite-time tomograms as in ( 1%T|) . For all different 0's 
the U(6), of which B (8) is the self-adjoint generator, are unitarily equivalent 
operators, hence all the tomograms share the same information. 

First we select a subset X n in such a way that the corresponding family 

1 ^e x W r ^ s orthogonal and normalized, 

< 4%AL >= <w (ss) 

This is possible by taking the sequence 

X n = X + — sin^ n e Z (84) 

where X is freely chosen (in general we take X = 0). We then consider the 
projections of the signal f(t) 

4„ (/)=</, 41 > (85) 

Denoising consists in eliminating the c x (/) such that 

K(/)| a <e ( 86 ) 
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for some threshold e. This power selective denoising is more robust than, 
for example, frequency filtering which may also eliminate important signal 
information. 

The component separation technique is based on the search for an inter- 
mediate value of 9 where a good compromise might be found between time 
localization and frequency information. This is achieved by selecting subsets 
Tk of the X n and reconstructing partial signals (fc-components) by restricting 
the sum to 

/*(*) = E c x n (f)^,x n (t) (87) 

n£T k 

for each k. 

As an example consider the following signal 

y{t)=y 1 (t)+y 2 {t)+y 3 (t) + b{t) (88) 

j/j (£) = exp (z25*) ,te[Q, 20] 

y 2 (t) = exp (i75t) , t e [0, 5] 

y 3 (t) = exp (z75t) , t E [10, 20] (89) 

Separation is impossible both at the time (9 = 0) and the frequency 
(9 = | ) axis. However, at some intermediate 9 value one obtains distinct 
probability peaks (Figj9]), which after the projections ( 187)) allows an accurate 
separation of the signal components (FigsHO] and ITT)) 

Component decomposition of more complex signals (nonlinear chirps over- 
lapping in both the time and the frequency domains and experimental re- 
flectometry signals) has been successfully carried out by this technique [26] 

E7J. 



8.3 Tomograms and turbulent velocity fluctuations 

Here we report briefly on an analysis by the tomographic technique of a ve- 
locity fluctuation signal of a turbulent flow in a wind tunnel. It illustrates 
the fact that the choice of the pair of non-commuting operators in tomogram, 
should be adapted to the signal under study. As before we use finite-time 
tomograms in the interval (to^o + T). For the finite-time (time- frequency) 
tomogram M 1; the normalization and a set of A n 's leading to an orthonor- 
malized set of eigenstates has already been written in (18T)) - (182l . 
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Figure 9: The tomogram at 9 — 7r/5 used for the component separation 

For future reference we include here the corresponding sets of orthonor- 
malized eigenstates for the finite-time time-scale tomogram M2(/i, u, X) (Eqj30 
and for the finite-time time-conformal tomogram M^fi, u,X) (Eq l32l : 



M 2 {9,X) 



t +T 2 

f*(t)4%(t)dt 



to 



4% (t) 



i 



i 



Vlog|t + T|-log|t | V\i 

2nn 



expi 



= \<f,4 2) > 
cos 9 X 



* 



and 



X n — Xo + 
M 4 (9,X) = 

4% (t) = 



log|* + T| -log |t 1 

to+T 



sin 9 sin 9 

sin 9 n G Z 



(90) 
log |*| } (91) 



,( 4 ) 



rmr x (t)dt 



to 



\<f,4 i} > 



t (t + T) 1 .(cos9 X 

expt ( — — - log |*| -f 



T t 



sin# 



*sin# 



(92) 

(93) 
(94) 
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Figure 10: The yi component 



A, 



, t (t + T) n . 
X H 27m sin # 



neZ 



(95) 



Analyzing the turbulent velocity fluctuations signal with these tomo- 
grams, one notices that except for some features on the frequency axis corre- 
sponding to some dominating frequencies, no interesting structures are put 
into evidence when one use the time-frequency tomogram. The situation is 
more interesting for the time-scale tomogram M2 (6,X). In Fig{12]we show 
a contour plot for M 2 (9, X) corresponding to a section of 1000 data points. 
For intermediate regions of 9 one notices, a strong concentration of energy in 
a few regions. This is put into evidence by a cut at 9 = 1.26 (Fig JT3|) . Pro- 
jecting out the signal corresponding to these regions with the corresponding 
rex W s a ^ thi s Q-i one sees that although the signal has many complex fea- 
tures most of the energy is concentrated in fairly regular structures. FigJT 
shows the structure 77 (£) corresponding to the second peak in FigJTSl 
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Figure 11: The 1/3 component 

9 Conclusions 

Tomograms provide a two-variable characterization of signals which, due to 
its rigorous probabilistic interpretation, is robust and free of artifacts and 
ambiguities. For each particular signal that one wants to analyse the choice 
of the appropriate tomogram depends not only on the signal but also on 
the features that we might want to identity or emphasize. So far we have 
explored component separation, denoising and identification of small signal 
in noise, but other features may also benefit from the robust probabilistic 
of the tomographic analysis. This was our main motivation to include here 
a long list of many different operator choice leading to different classes of 
tomograms. 

The description of the tomograms as operator symbols, with the cor- 
responding quantizers and dequantizers, not only provides an alternative 
formulation but may also be used to extend the algebraic signal processing 
formalism to a wider nonlinear context. 
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Figure 12: Contour plot of the tomogram for a velocity fluctuations signal 
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